Three-potential formalism for the three-body scattering problem with attractive 

Coulomb interactions 
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A three-body scattering process in the presence of Coulomb interaction can be decomposed formaUy 
into a two-body single channel, a two-body multichannel and a genuine three-body scattering. The 
corresponding integral equations are coupled Lippmann-Schwinger and Faddeev-Merkuriev integral 
equations. We solve them by applying the Coulomb-Sturmian separable expansion method. We 
present elastic scattering and reaction cross sections of the e'^ + H system both below and above 
the H{n — 2) threshold. We found excellent agreements with previous calculations in most cases. 
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The three-body Coulomb scattering problem is one 
of the most challenging long-standing problems of non- 
relativistic quantum mechanics. The source of the dif- 
ficulties is related to the long-range character of the 
Coulomb potential. In the standard scattering theory 
it is supposed that the particles move freely asymptoti- 
cally. That is not the case if Coulombic interactions are 
involved. As a result the fundamental equations of the 
three-body problems, the Faddeev-equations, become ill- 
behaved if they are applied for Coulomb potentials in a 
straightforward manner. 

The first, and formally exact, approach was proposed 
by Noble His formulation was designed for solv- 

ing the nuclear three-body Coulomb problem, where all 
Coulomb interactions are repulsive. The interactions 
were split into short-range and long-range Coulomb-like 
parts and the long-range parts were formally included in 
the "free" Green's operator. Therefore the corresponding 
Faddeev-Noble equations become mathematically well- 
behaved and in the absence of Coulomb interaction they 
fall back to the standard equations. However, the asso- 
ciated Green's operator is not known. This formalism, 
as presented at that time, was not suitable for practical 
calculations. 

In Noble's approach the separation of the Coulomb- like 
potential into short-range and long-range parts were car- 
ried out in the two-body configuration space. Merkuriev 
extended the idea of Noble by performing the splitting 
in the three-body configuration space. This was a cru- 
cial development since it made possible to treat attrac- 
tive Coulomb interactions on an equal footing as repul- 
sive ones. This theory has been developed using in- 
tegral equations with connected (compact) kernels and 
transformed into configuration-space differential equa- 
tions with asymptotic boundary conditions In prac- 
tical calculations, so far only the latter version of the 
theory has been considered. The primary reason is that 
the more complicated structure of the Green's operators 
in the kernels of the Faddeev-Merkuriev integral equa- 



tions has not yet allowed any direct solution. However, 
use of integral equations is a very appealing approach 
since no boundary conditions are required. 

Recently, one of us has developed a novel method for 
treating the three-body problem with repulsive Coulomb 
interactions in three-potential picture ^J. In this ap- 
proach a three-body Coulomb scattering process can be 
decomposed formally into a two-body single channel, a 
two-body multichannel and a genuine three-body scat- 
tering. The corresponding integral equations are coupled 
Lippmann-Schwinger and Faddeev-Noble integral equa- 
tions, which were solved by using the Coulomb-Sturmian 
separable expansion method. The approach was tested 
first for bound-state problems Q with repulsive Coulomb 
plus nuclear potential. Then it was extended to calculate 
p — d scattering at energies below the breakup threshold 
1^ and more recently we have used the method to cal- 
culate resonances of three- a systems jsj. Also atomic 
bound-state problems with attractive Coulomb interac- 
tions have been considered ||] . These calculations showed 
an excellent agreement with the results of other well es- 
tablished methods. The efficiency and the accuracy of 
the method was demonstrated. 

The aim of this paper is to generalize this method 
for solving the three-body Coulomb problem with re- 
pulsive and attractive Coulomb interactions. We com- 
bine the concept of three-potential formalism with 
the Merkuriev's splitting of the interactions and solve 
the resulting set of Lippmann-Schwinger and Faddeev- 
Merkuriev integral equations by applying the Coulomb- 
Sturmian separable expansion method. In this paper 
we restrict ourselves to energies below the three-body 
breakup threshold. 
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I. INTEGRAL EQUATIONS OF THE 
THREE-POTENTIAL PICTURE 

We consider a three-body system with Hamiltonian 

F = i/" + «f + + (1) 

where H'^ is the three-body kinetic energy operator and 
denotes the Coulomb-Uke interaction in subsystem 
a. The potential may have repulsive or attractive 
Coulomb tail and any short-range component. We use 
the usual configuration-space Jacobi coordinates Xa and 
ya] Xa is the coordinate between the pair (/?, 7) and ya 
is the coordinate between the particle a and the center 
of mass of the pair 7). Thus the potential v^, the 
interaction of the pair (/3,7), appears as v^{xa)- We 
also use the notation X = {xa,ya} G R^- 

A. Merkuriev's cut of the Coulomb potential 

The Hamiltonian (Q) is defined in the three-body 
Hilbert space. The two-body potential operators are for- 
mally embedded in the three-body Hilbert space 

v'' ^v^^ix^y. (2) 

Merkuriev introduced a separation of the three-body con- 
figuration space into different asymptotic regions. The 
two-body asymptotic region Qa is defined as a part of 
the three-body configuration space where the conditions 

K|<:rO(l + |y„|/y°)l/^ (3) 

with x^jVa > and v > 2, are satisfied. He proposed to 
split the Coulomb interaction in the three-body configu- 
ration space into short-range and long-range terms 

v^^vi^^+v(^\ (4) 

where the superscripts s and / indicate the short- and 
long-range attributes, respectively. The splitting is car- 
ried out with the help of a splitting function 

v<''Hx,y)^v^{x)ax,y), (5) 
v^'Hx,y)^v''{x)[l-ax,y)]. (6) 

The function ( is defined such that 

ax,y) otherwise. 

In practice, in the configuration-space differential equa- 
tion approaches, usually the functional form 

ax,y)^2/{l+[ix/x"r/il + y/y')]}, (8) 

was used. 



The long-range Hamiltonian is defined as 

i7«=i70 + t-«+4')+z;W, (9) 

and its resolvent operator is 

G«(z) = (10) 

Then, the three-body Hamiltonian takes the form 

H ^ H^'^ + vi'^ + v],'^ + v^^'l (11) 

In the conventional Faddeev theory the wave function 
components are defined by 

\^a)^{z-H°)-^Va\^), (12) 

where is a short-range potential and is the Fad- 
deev component of the total wave function j^I'). While 
the total wave function j^*), in general, has three dif- 
ferent kind of two-body asymptotic channels, \tpa) pos- 
sesses only a-type two-body asymptotic channel. The 
other channels are suppressed by the short-range poten- 
tial Va ■ This procedure is called asymptotic filtering and 
it guarantees the asymptotic orthogonality of the Fad- 
deev components [Q. 

The aim of the Merkuriev procedure was to formally 
obtain a three-body Hamiltonian with short-range poten- 
tials w'-'*^ and long-range Hamiltonian i?*^'^ in order that 
we can repeat the procedure of the conventional Faddeev 
theory. The total wave function \'^) is split into three 
components, 

I*) = Ha) + 1^/3) + 1^7). (13) 

with components defined by 

IV^a) =GW«i^)|*). (14) 

This procedure is an example of asymptotic filtering. The 
short-range potential Va ^ acting on j^") suppresses the 
possible f3 and 7 asymptotic two-body channels, pro- 
vided G*^'^ itself does not introduce any new two-body 
asymptotic channels. With the Merkuriev splitting this 
is avoided because iJ^'-* does not have two-body asymp- 
totic channels even if some of the long-range potentials 
have attractive Coulomb tail. In the attractive case w*-'^ 
appears as a valley along the y — x'^ parabola-like curve 
with Coulomb-like asymptotic behavior in x at any finite 
y. (See Figs, ^and ^for the short- and long-range parts, 
respectively). However, as ?/ — > 00 the depth of the valley 
goes to zero, consequently the two-body bound states are 
pushed up, and finally the system does not have any two- 
body asymptotic channels. We note that the Merkuriev 
formalism contains the Noble's in the limit y^ — > 00. 

B. The three-potential picture 

In Ref. 1^ the three body scattering problem with 
repulsive Coulomb interactions were considered in the 
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three-potential picture. In this picture the scattering 
process can be decomposed formally into three consec- 
utive scattering processes: a two-body single channel, a 
two-body multichannel and a genuine three-body scatter- 
ing. This formalism also provides the integral equations 
and the method of constructing the S'-matrix. Below we 
adapt this formalism to attractive Coulomb interactions 
along the Merkuriev approach. 

The asymptotic Hamiltonian is defined as 

H^=H°+v^, (15) 

and the asymptotic states are the eigenstates of Ha 

(16) 

where {xaya\^a) — {ya\Xa){xa\<t'a) is a product of a 
scattering state in coordinate ya and a bound state in 
the two-body subsystem Xa- 

We define the two asymptotic long-range Hamiltonians 

as 

h2^^H' + v^ + v^^^+v\'^ (17) 

and 

Ha=H° + v^ + u^S', (18) 

where Ua"* is an auxiliary potential in coordinate ?/q,, and 
it is required to have the asymptotic form 

-Z„(Z^ + Z^)/2/, (19) 

as ya — oo. In fact, is an effective Coulomb-like 
interaction between the center of mass of the subsystem 
a (with charge Zfs + Z-y) and the third particle (with 
charge Za). We introduced this potential in order that 
we compensate the long range Coulomb tail of vj^^ + v^"^ 
in Qa- 

Let us introduce the resolvent operators: 

G{z)^iz-Hy\ (20) 

GW(^) = (z-i?('))-\ (21) 

Ga{z) ^ (z ~ Ha)-\ (22) 

The operator Gq"* is the long-range channel Green's oper- 
ator and Ga is the channel distorted long-range Green's 
operator. These operators are connected via the follow- 
ing resolvent relations: 

G{z) = Gli\z) + gW(z)1^«G(z), (23) 
GW(z) = Ga{z) + G„(z)(7"Gi')(z), (24) 



where V" = u^"^ -t- u^^"^ and = vf + ^ 

The scattering state, which evolves from the asymp- 
totic state \<^a) under the influence oi H, is given as 

l^-W) = limieG(£;„±i£)|$„}. (25) 

Similarly, we can define the following auxiliary scattering 
states 

= limieG«(ii;±i£)|$„) (26) 

and 

= lim ieG,(£; ± i£)|$a), (27) 

which describe scattering processes due to Hamiltonians 
Ha'' and Ha, respectively. 

The S-matrix elements of scattering processes are ob- 
tained from the resolvent of the total Hamiltonian by the 
reduction technique |^ 

SpJ^a^ = ^hm Ihn ieei(^«-^")*($0,|G(£;c., + ie)\^a^)■ 

(28) 

The subscript i and j denotes the i-ih and j-th eigen- 
states of the corresponding subsystems, respectively. If 
we substitute (^3|) into (|2^) we get the following two 
terms: 

limieei(^«-^")*($/3,|GL'ni?m+ie)l'i>m) 

(29) 

4?- = /^m limi£ei(^«-^")*($^,|GW(£;„ -f i£) 

V''G{Ea^+ie)\^^a^). (30) 

Substituting Eq. ( p^ ) into (|29|), the first term yields two 
more terms 

4Va,= 1™ limieei(^«-^"'*($0,|G„(£;„ + i£)|$„) 

(31) 

= lim lim ieei(^«-^=')*(<D^,IG.(i?.. -t-ie) 

U''G'^l\Ea^+ie)\^c.^). (32) 

Using of the properties of the resolvent operators the 
limits can be performed and we arrive at the following, 
physically plausible, result. The first term, S^^j^^, is the 
S-matrix of a two-body single channel scattering on the 
potential 

= 5pa5,.S{u<{^). (33) 

If Ma'' is a pure Coulomb interaction S{u"a) falls back 
to the S-matrix of the Rutherford scattering, if u^a is 
identically zero s'g^ equals to unity. The second term. 
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^(2) 



describes a two-body multichannel scattering on 



the potential 

^So. = "27ri<50„<5(£;^, - Eo.^){^-^ |C/"l<i'S^''^ (34) 

The third term gives account of the complete three-body 
dynamics 



s: 



(3) 

Pj,ai 



2niSiEp,~E^,){<^>fl-^\V'^\¥^f). (35) 



C. Lippmann-Schwinger integral equation for 1$^^) 

Starting from the definition of |<&a^), Eq. (p6|), by uti- 
lizing the resolvent relation ( |2^ ) and the definition (27) 
we easily derive a Lippmann-Schwinger equation 



= |<fi±)) +Ga{E± ie)C/"|$«W), (36) 
where are given by 

I^L^^) = lxL^^)l0a>- (37) 

The state Ixi^"*} is a scattering state in the Coulomb-like 
potential Ua\ya)- 



Merkuriev showed that after a certain number of iter- 
ations these equations were reduced to Fredholm inte- 
gral equations of the second kind with compact kernels 
for all energies, including energies below (E < 0) and 
above {E > 0) the three-body breakup threshold 
Thus all the nice properties of the original Faddeev equa- 
tions established for short-range interactions remain valid 
also for the case of Coulomb-like potentials. We note 
that the triad of Lippmann-Schwinger equations and the 
set of Faddeev equations describe the same physics, the 
equations have identical spectra and in fact, the Faddeev 
equations are the adjoint representations of the triad of 
Lippmann-Schwinger equations pO| ]. 

Utilizing the properties of the Faddeev components the 
matrix elements in ( |35| ) can be rewritten in a form better 
suited for numerical calculations 



(0(- 



(')(-)!„, («)|„;,(+) 



/3j 



(44) 



Summarizing, in the three-potential formalism, start- 
ing from l^i^-*), by solving a Lippmann-Schwinger equa- 



tion, we determine |$a'*''^' 



Then from by 



solving the set of Faddeev-Merkuriev equations, we de- 
termine the components |V'a^^)- Finally using Eqs. ( |3^ ) 
and (H4) we construct the iS-matrix. 



D. Faddeev-Merkuriev integral equations for the 
wave function components 

The integral equations for the wave function I^'q^'') are 
arrived at by combining the resolvent relation ( |2^ ) and 
Eq. (^^. In this case however we have three resolvent 
relations and therefore we obtain a triad of Lippmann- 
Schwinger equations 

l^-W) ^ |$W(±)^ + GW(^±iO)y"|*W) (38) 

G^''(£;±iO)l/''|*i±)) (39) 

Gij\E±iO)V^\-^t'>). (40) 

Although these three equations together provide unique 
solutions their kernels are not connected therefore 
they cannot be solved by iterations. The way out of 
the problem is to use the Faddeev decomposition which 
leads to equations with connected kernels, thus they are 
effectively Fredholm- type integral equations. 

Multiplying each elements of the triad from left by 
G'^''^Va^ and utilizing ( p^ ) we get the set of Faddeev- 
Merkuriev integral equations for the components 

IV^W) = |$»(±)> + G^^He ±iO)vin4^^) + 1^7^^)] (41) 
1^^?^)= Gy(i?±iO)z;(^)[|V^W)-)-|^i±))] (42) 

|V^^))= Gi'\E±iO)vi;^[\^i^^) + \A^^)]. (43) 



II. COULOMB-STURMIAN SEPARABLE 

EXPANSION APPROACH TO THE 
THREE-BODY INTEGRAL EQUATIONS 

In order to solve operator equations in quantum me- 
chanics one needs a suitable representation for the oper- 
ators. For solving integral equations it is especially ad- 
vantageous if one uses a representation where the Green's 
operator is simple. For the two-body Coulomb Green's 
operator there exists a Hilbert-space basis in which its 
representation is very simple. This is the Coulomb- 
Sturmian (CS) basis. In this representation-space the 
Coulomb Green's operator can be given by simple and 
well-computable analytic functions jll] . This basis forms 
a countable set. If we represent the interaction term on 
a finite subset of the basis it looks like a kind of sep- 
arable expansion of the potential, and so the integral 
equation becomes a set of algebraic equations which can 
then be solved without any further approximation. The 
completeness of the basis ensures the convergence of the 
method. 

This approximation scheme has been thoroughly tested 
in two-body calculations. Bound- and resonant-state cal- 
culations were presented first ||Tl[| . Then the method 
was extended to scattering states Q. Since only the 
asymptotically irrelevant short-range interaction is ap- 
proximated, the correct Coulomb asymptotic is guaran- 
teed ]l3[. A recent account of this method is presented 
in Rcf. ]14|] . The method also proved to be very efficient 
in solving three-body Faddeev-Noble integral equations 
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for bound- |Q] and scattering-state jsj problems with re- 
pulsive Coulomb interactions. 

In subsection A we define the basis states in two- and 
three-particle Hilbert space. In subsection B we review 
some of the most important formulae of the two-body 
problem. In subsections C and D we describe the calcu- 
lation of the S'-matrix and the solution of the Faddeev- 
Merkuriev integral equations. We follow the line pre- 
sented in Ref. S. 



and consider the inhomogeneous Lippmann-Schwinger 
equation for the scattering state [ipi) in some partial wave 



I 



lv^^) = l0p)+.9F(i?)«PV^)- 



(50) 



Here \4>f') is the regular Coulomb function, gf[E) is the 
two-body Coulomb Green's operator 



C\-l 



(51) 



A. Basis states 

The Coulomb-Sturmian functions in some angular 
momentum state I are defined as 

1/2 

{r\nl) 



(n + 2Z + l)! 



(26r)'+i exp(-6r)L2'+i(26r), 

(45) 

n = 0, 1, 2, . . . . Here, L represents the Laguerre polyno- 
mials and 6 is a fixed parameter. In an angular momen- 
tum subspace they form a complete set 



N 



1 = lim > \nl){nl\ = lim Iw, 

JV^oo ^ N^oc 



(46) 



where \nl) in configuration-space representation reads 
{r\nl) — {r\nl)/r. 

The three-body Hilbert space is a direct sum of two- 
body Hilbert spaces. Thus, the appropriate basis in an- 
gular momentum representation should be defined as a 
direct product 

\nvl\)a^\nl)c,®\v\)c,, (n,iy = 0,1,2,...), (47) 



with the CS states of Eq. ( |45| ) . Here I and A denote the 
angular momenta associated with Jacobi coordinates x 
and y, respectively. In our three-body Hilbert space basis 
we take bipolar harmonics in the angular variables and 
CS functions in the radial coordinates. The complete- 
ness relation takes the form (with angular momentum 
summation implicitly included) 



1 = fim 

N- 



N 

,E 

n,f — 



(48) 



where {xaya\'nh'l^)a = {xaVal nvlX)a/[xaya)- It should 
be noted that in the three-particle Hilbert space we can 
introduce three equivalent basis sets which belong to frag- 
mentation a, /3 and 7. 



with the free Hamiltonian h^. We make the following 
approximation on Eq. (pC 



Wi) = \v?)+9?{E)lNv\'hN\4^i) 



(52) 



is) 

i.e. we approximate the short-range potential by a 
separable form 



N 



= lim lArvf^'lTv 



n,n'— 



where the matrix 



is) 



(53) 



(54) 



These matrix elements can always be calculated (numer- 
ically) for any reasonable short-range potential. In prac- 
tice we use Gauss-Laguerre quadrature, which is well- 
suited to the CS basis. _ 

Multiplied with the CS states {nl\ from the left, 
Eq. (52) turns into a linear system of equations for the 
wave- function coefficients ip = {nl\ipi) 



(55) 



where the underlined quantities are matrices with the 
following elements 



and 



V? = {nl\^?) (56) 



gf ^iE)^{nl\gf{E)\n'l). (57) 



1. The matrix elements (nl\gf (z)\n' I) 



B. Coulomb-Sturmian separable expansion in 
two-body scattering problems 

Let us study a two-body case of short-range plus 
Coulomb-like interactions 



vi 



(49) 



The key point in the whole procedure is the exact 
and analytic calculation of the CS matrix elements of 
the Coulomb Green's operator and of the overlap of the 
Coulomb and CS functions. For the Green's matrix we 
have developed two independent, analytic approaches. 
Both are based on the observation that the Coulomb 
Hamiltonian possesses an infinite symmetric tridiagonal 
(Jacobi) matrix structure on CS basis. 
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Let us consider the radial Coulomb Hamiltonian 

^2 / j2 



' 2m V dr^ 



r 



(58) 



where m, Z and Z stands for the mass, angular mo- 
mentum and charge, respectively. The matrix J^„/ = 
{n\{z — h'j^)\n') possesses a Jacobi structure. 



J^^ = 2{n + l + l){k^-b') 



- Z 



and 



Amh ' 



(59) 



(60) 



where k = {2mz /'ti?Y^'^ is the wave number. The main 
result of Ref. flq] is that for Jacobi matrix systems the 
A^'th leading submatrix g^^^ of the infinite Green's ma- 
trix can be determined by the elements of the Jacobi 

matrix 



Replacing the tail of the continued fraction by its explicit 
analytical form w±^ we can speed up the convergence 
and, more importantly turn a non-convergent continued 
fraction into a convergent one . Analytic continuation 
is achieved by using w± instead of the non-converging 
tail. In Ref. 161 it was shown that w+ provides an ana- 
lytic continuation of the Green's matrix to the physical, 
while i/7_ to the unphysical Riemann-sheet. This way 
Eq. ( |6^ ) together with (|6l|) provides the CS basis repre- 
sentation of the Coulomb Green's operator on the whole 
complex k plane. We note here that with the choice of 
Z = the Coulomb Hamiltonian ( |58| ) reduces to the ki- 
netic energy operator and our formulas provide the CS 
basis representation of the Green's operator of the free 
particle as well. We emphasize that this procedure does 
not truncate the Coulomb Hamiltonian, because all the 
higher J„„/ matrix elements are implicitly contained in 
the continued fraction. 

We note that g'~' has already been calculated before 
]lT[ . From the J-matrix structure a three-term recursion 
relation follows for the matrix elements . This recur- 

— Tin' 

sion relation is solvable if the first element q'-^ is known. 

—00 
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^'■"I^-' = [J_nn' + SnN ^n'N dSiN+i ^] ^ i (61) It is given in a closed analytic form 



where C is a continued fraction 

UN 



C=-- 



dN 



UN+l 



dN+1 



UN+2 



dN+2 + 



(62) 



Amh 1 



^00 



{b-ikfl + irj+1 
'X2F1 -l + bil;l + iri + 2, 



b + ik 
b-ik 



, (69) 



with coefficients 



ii-n,n-l/ 



LL.n,nl —1 



■n,n+l 



(63) 



In Ref. |T|] it was shown that although the continued 
fraction C is convergent only on the upper-half k plane it 
can be continued analytically to the whole k plane. This 
is because the and d„ coefficients satisfy the limit 
properties 



u — lim Un = — 1 

n — >oo 



d = lim dn = 2{k^ ~ b^)/ik^ + b^ 

ri — >oo 

Then the continued fraction appears as 

UN 



C 



dN 



UN+l 



(64) 
(65) 

(66) 



"TF 



Therefore the tail w of C satisfies the implicit relation 

u 



d + w 

which is solved by 

w± = {b±ikf/{b^ + k^) 



(67) 



(68) 



where rj = Zm/ {Ti k) is the Coulomb parameter and 2-P1 
is the hypergeometric function. For those cases where the 
first or the second index of 2^1 is equal to unity, there 
exists a continued fraction representation, which is very 
efficient in practical calculations. It was shown that the 
two methods lead to numerically identical results for all 
energies and our numerical continued fraction represen- 
tation possesses all the analytic properties of g'~^ . The 
exact analytic knowledge of g'-^ allows us to calculate the 
matrix elements of the full Green's operator in the whole 
complex plane 



(70) 



The overlap vector of CS and the Coulomb functions 
{nl\(pf) is known analytically It can be calculated 
by a three-term recursion, derived from the J-matrix, us- 
ing the starting value 



(OZ|(pp) = cxp(277 arctan(fc/6))4 



2TTrj 



exp(27rr/) — 1 



2k/b 



i+i I 



n 



1/2 



(71) 
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C. Calculation of the three-body S-matrix 



The aim of any scattering calculation is to determine 
the S-matrix elements. In our case we need to calcu- 
late the terms (||), (||) and (|4|) of the three-potential 
picture. 

The term S^g^^^ is trivial because it is just the two- 



matrix elements {^[^j''\U°'\^ 
a-) 



(*^Q7 I contains a two-body bound-state wave function 



in coordinate Xa this matrix element is confined to fla, 
where C/" is of short-range type. Therefore a separable 
approximation is justified 



We conclude that to calculate the S-matrix of the 
three-potential formulae we need the CS matrix elements 
( [7^ ) and (^), which can always be evaluated numeri- 
cally by using the transformation of Jacobi coordinates 
jlSf . In addition we need the CS wave function com- 
ponents a{ni^lX\^^J'i''), a{ni^lM^aP) ^nd a{ni^lX\ip^o^^)- 
We determine them in the following section by solv- 
ing Lippmann-Schwinger and Faddeev-Merkuriev inte- 
gral equations. 

It should be noted that the approximations (173) and 
} are needed. Since (76) used in calculating the matrix elements ( [7q ) and 

(78) become equalities as N goes to infinity. In practical 
calculations we increase N until we observe numerical 
convergence in scattering observables. 



body S-matrix of the Coulomb-like potential Ua^. 

To calculate the second term, S'i^-^^- of Eq. (|3^), the 



(-) 



), (72) 



D. Solution of the three-body integral equations 



i.e, in this matrix element, we can approximate U" by a 
separable form 



In the set of Faddeev-Merkuriev equations (^l|-^3|) we 
make the approximation of ( [76| ) 



iJ" = lim l^i7"l^ 



1a jja-i a 



N 

E 

n.u.n' m'—O 



(73) 



IV-a) = + G«[l^«('')4|#) + l^^i^^)!^!^,)] (79) 
Gg)[<4^h;|V^^) + 44^)l^|V'a)] (80) 
1^,}= GW[1>(^)i^|^„) + 1>(^)<|V'/5)]. (81) 



where 



■nulX.n' u 



I'X' 



a{nvlX\U°'\nvl'X) 



(74) 



The matrix element appears as 



N 



\ OLJ I \ Oil 



(0(^ 



«(+)\ 
(75) 



In calculating the third term, 5*^^'^^, of (|44|), we have 
matrix elements of the type ($^7'' Iwi^"* IV'^t'')- Here we 



(s) 

can again approximate the short-range potential Va in Qj^d 
the three-body Hilbert space by a separable form 



N 



nni ijv^^Q i-N 



N 



^ \nvl\)^^^lsi{n'v'l'\'\ 

n,i',n' ,i''—0 



(76) 



where 



Multiplied by the CS states c((nz//A|, ^(nz^ZAj and^(ni^/A|, 
respectively, from the left the set of integral equations 
turn into a linear system of algebraic equations for 



the coefhcients of the Faddeev components -0 



(0^-l<hm 



(82) 



with 



G(') 



{niylX\G'^!^^\n'iy'l'X')o,, (83) 



$(0 



.a|a>(')>. 



(84) 



Notice that the matrix elements of the Green's operator 
are needed only between the same partition a whereas 
the matrix elements of the potentials occur only between 
different partitions a and /3. 



a{ni^lX\vi'^\n'iy'l'X')fj. (77) 



In ( |7^ ) the ket and bra states belong to different frag- 
mentations depending on the neighbors of the potential 
operators in the matrix elements. Finally, the matrix 
elements take the form 



N 



IniylX}^^/^^ p{n'iy'l'X'\^'^^ 
(78) 



/,(+)\ 



1. The matrix elements a{nul\\Ga^ \n' u'l' \')a and 
c«(n1WA|$L''> 

Unfortunately neither the matrix elements ( |8^ nor the 
overlaps ( p^ are known. The appropriate Lippmann- 
Schwinger equation for Ga^ was proposed by Merkuriev 



G<^\z) = G-(2) + G-(z)yrG('H^), (85) 
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where G^* and V^* are the asymptotic channel Green's 
operator and potential, respectively. A similar equation 
IS valid for 



i$W) = i$-> + G-(z)Fri$«). 



(86) 



Both gL'' and are genuine three-body quantities. 

One may wonder why a single Lippmann-Schwinger equa- 
tion suffices. The Hamiltonian Ha'' has a peculiar prop- 
erty - it has only a- type two-body asymptotic channels. 
For such systems a single Lippmann-Schwinger equation 
provides a unique solution |19|. 

The objects G"'', V^^ and <i>^'' are very complicated. 
Their leading order terms were constructed in configura- 
tions space in the different asymptotic regions. The po- 
tential V^^, as |X| ^ oo, decays faster than the Coulomb 
potential in all directions of the three-body configuration 
space: ~ 0{\X\-^-''), e > @. Therefore we may 
express the solutions of Eqs. (p^ and (p^) formally as 



and 



(87) 



respectively, where 

Go!,,,, -a(ni^/A|G^^|nVrA')„, 



(89) 



and 



(91) 



Here, G^", V^^ and appear between finite number of 
of square-integrable CS states, which confine the domain 
of integration to ^la- In this region, however, G^^ coin- 
cides with G„, V^' with U" and with §. Finally 
we have 



(G«)-i = (GJ~i-C/«, 



(92) 



where 



and 



G„ 



a{nh'lX\Ga\n'iy'l'\')a (93) 



U-Lix.n'u'i'X' = a{nvl\\U°'\n'v'l'X'), 
And in a similar way 

[(GJ-i = (GJ-ii„, 

where 



a{nvl\\<^a)- 



(94) 
(95) 
(96) 



We note that from Eq. (|9^) follows that the left side 
of Eq. ( |95| ) is just the inhomogeneous term of Eq. (82). 
Both Eqs. (|9^) and (|8^) are solved with the same inho- 
mogeneous term. 



2. The matrix elements a{'nvl\\Ga\n' v'l' \')a and 

The three-particle free Hamiltonian can be written as 
a sum of two-particle free Hamiltonians 



(97) 



Then the Hamiltonian Ha of Eq. (|l8| ) appears as a sum 
of two Hamiltonians acting on different coordinates 



Ha — 



(98) 



with — h^xa + '^ai^a) and hy 
which, of course, commute. The state l^^), which is 
an eigenstate of Ha , is a product of a two-body bound- 
state wave function in coordinate Xa and a two-body 
scattering-state wave function in coordinate j/q. Their 
CS representations arc known from the two-particle case 
described before. _ 

The matrix elements of Ga can be determined by mak- 
ing use of the convolution theorem 

Ga{z) = [z - hx^ - hy^y^ 

" 2^ / "^^'^^ " '^y-^^'- 

The contour G should encircle, in positive direction, the 
spectrum of hy^ without penetrating into the spectrum 
of hx^ ■ 

The convolution theorem follows from a more general 
formula. A function of a self adjoint operator h is defined 
as 



fih) 



1 
27ri 



dzf{z){z — h) 



(100) 



where G is a contour around the spectrum of h and / 
should be analytic on the region encircled by G. 

In the following we suppose that u^^^ either vanishes 
or is a repulsive Coulomb-like potential. This assump- 
tion is not necessary but it greatly simplifies the analysis 
below. Numerical examples show that there are a great 
many physical three-body systems where this condition 
is satisfied. This condition ensures that hy does not have 
bound states. 

To examine the analytic structure of the integrand (|99| ) 
let us shift the spectrum of gx^ by taking z ^ E + le 
with positive e. In doing so, the two spectra become 
well separated and the spectrum of gy^ can be encircled. 
The contour G is deformed analytically in such a way 
that the upper part descends to the unphysical Riemann 
sheet of gy^ , while the lower part of G can be detoured 
away from the cut [see Fig. ^ . The contour still encircles 
the branch cut singularity of gy^ , but in the e — > limit 
avoids the singularities of gx^- Thus, the mathemati- 
cal conditions for the contour integral representation of 
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Gq,(z) in Eq. ( P9[ ) is met. The matrix elements can 
be cast in the form 



1 



27ri 



Gaiz) = TT^ <P dz' g {z - z') g (z'), 



(101) 



where the corresponding CS matrix elements of the two- 
body Green's operators in the integrand are known ana- 
lytically for all complex energies. 



III. TEST OF THE METHOD 

We demonstrate the power of this new method by cal- 
culating elastic phase shifts oi + H scattering below 
the Ps{n = 1) threshold and cross sections of the e+ + 
elastic scattering as well as j<+ + Ps reaction channels up 
to the Ps{n = 2) threshold. In all examples we have to- 
tal angular momentum L — Q and we have taken angular 
momentum channels up to Z = 10. We use atomic units. 

Let us numerate the particles e+, p and e~, with 
masses me± = Ime and nip = 1836.1527me, by 1, 2 and 
3, respectively. In the channel 3 there are no two-body 
asymptotic channels since the particles e+ and p do not 



form bound states. Therefore, we can take v 



and 



include the total Wg in the long range Hamiltonian 



H = ij(') + 



„(0 



(102) 
(103) 



In this case j^/ia) = and we have the set of two- 
component Faddeev-Merkuriev equations 



1^2) 



IV'2). 



(104) 
(105) 



The parameters of the splitting function C, of Eq. (g) are 
rather arbitrary. The final converged results should be 
insensitive to their values; our numerical experiences con- 
firm this expectation. For the parameters of C we have 
taken v = 2.1, — Z and — 10, whereas for the 
parameters of CS functions we have taken h = 0.9. We 
have experienced that the rate of convergence is rather 
insensitive on the choice of h over a broad interval. 

First we examine the convergence of the results for 
cross sections at incident wave numbers fci = 0.71, 
fci = 0.75 and fci — 0.8, which correspond to scatter- 
ing states in the Ore gap. Table | shows the conver- 
gence of e+ + H— > + H elastic scattering (an) and 
e+ + H— > p+ + Ps positronium formation (CT12) cross 
sections (in ttuq) with respect to N, the number of CS 
functions in the expansion, and with respect to increas- 
ing the angular momentum channels in the bipolar ex- 
pansion. For comparison we provide the results of Ref. 
pO[ . We can see that very good accuracy is achieved even 
with relatively low N in the expansion. 

In Table we compare our converged results for phase 
shifts (in radians) below the Ps{n — 1) threshold to that 



of other methods. Ref. |2^ is the best variational calcu- 
lation. In Ref. the Schrodinger equation was solved 
by means of finite-element method. In Refs. ||2^ and |2^] 
the configurations space Faddeev-Merkuriev differential 
equations were solved using the bipolar harmonic expan- 
sion method and in total angular momentum represen- 
tation, respectively. We can report perfect agreements 
with previous calculations. 

In Table III we present partial cross sections in 
the H{n = 2) - Ps{n = 2) gap (threshold energies 
0.7496-0.8745 Ry). In Ref. (2|] the configurations space 
Faddeev-Merkuriev differential equations were solved us- 
ing the bipolar harmonic expansion in the angular vari- 
ables an quintic spline expansion in the radial coordi- 
nates. We can report fairly good agreements. 



IV. CONCLUSION 

We have extended the three-potential formalism for 
treating the three-body scattering problem with all kinds 
of Coulomb interactions including attractive ones. We 
adopted Merkuriev's approach and split the Coulomb po- 
tentials in the three-body configuration space into short- 
range and long-range terms. In this picture the three- 
body Coulomb scattering process can be decomposed into 
a single channel Coulomb scattering, a two-body multi- 
channel scattering on the intermediate-range polarization 
potential and a genuinely three-body scattering due to 
the short-range potentials. The formalism provides us a 
set of Lippmann-Schwinger and Faddeev-Merkuriev in- 
tegral equations. 

These integral equations are certainly too complicated 
for the most of the numerical methods available in the 
literature. The Coulomb-Sturmian separable expansion 
method can be successfully applied. It solves the three- 
body integral equations by expanding only the short- 
range terms in a separable form on Coulomb-Sturmian 
basis while treating the long-range terms in an exact 
manner via a proper integral representation of the three- 
body channel distorted Coulomb Green's operator. The 
use of the Coulomb-Sturmian basis is essential as it allows 
an exact analytic representation of the two-body Green's 
operator, and thus the contour integral for the channel 
distorted Coulomb Green's operator can be calculated. 
The method provides solutions which are asymptotically 
correct, at least in fla, which is sufficient if the scattering 
process starts from a two-body asymptotic state. Since 
the two-body Coulomb Green's operator is exactly calcu- 
lated all thresholds are automatically in the right location 
irrespective of the rank of the separable approximation. 
The method possesses good convergence properties and 
in practice it can be made arbitrarily accurate by em- 
ploying an increasing number of terms in the expansion. 
Certainly, there is plenty of room for improvement but 
we are convinced that this method can be a very power- 
ful tool for studying three-body systems with Coulomb 
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interactions. 
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TABLE L Convergence of e + H— > e + H elastic scat- 
tering ((Til) and e"*" -|- H— > p + Ps positronium formation 
((712) cross sections (in tvOo) with respect to A'^, the number of 
CS functions in the expansion, and with respect to increasing 
the angular momentum channels (Imax) in the bipolar basis. 





^max 6 


^max 8 


^max 10 1 


iV 


0"11 (Tl2 


0"11 (Tl2 


(Til (Tl2 



12 


0.02662 


13 


0.02608 


14 


0.02581 


15 


0.02562 


16 


0.02548 


17 


0.02541 


18 


0.02532 


19 


0.02528 


20 


0.02522 



0.00423 
0.00424 
0.00423 
0.00424 
0.00425 
0.00426 
0.00427 
0.00427 
0.00428 



ki ^ 0.71, Ref. pc|]: an = 0.025, ai2 = 0.0038 



ki = 0.75, Ref. 



04412 
04345 
04318 
04280 
04269 
04252 
04246 
04238 
04232 



0.00441 
0.00440 
0.00440 
0.00440 
0.00440 
0.00441 
0.00442 
0.00442 
0.00442 



0.02664 0.00397 0.02665 0, 

0.02609 0.00398 0.02610 0, 

0.02582 0.00398 0.02583 0, 

0.02561 0.00398 0.02562 0, 

0.02546 0.00400 0.02547 0, 

0.02539 0.00401 0.02539 0, 

0.02529 0.00401 0.02530 0, 

0.02524 0.00402 0.02525 0, 

0.02517 0.00403 1 0.02518 0. 
pc|]: an = 0.044, (T12 = 0.0043 



00393 
00394 
00394 
00395 
00396 
00397 
00398 
00398 
00399 



ki = 0.80, Ref. 



06572 
06573 
06518 
06488 
06457 
06440 
06427 
06418 
06412 



0.00475 
0.00481 
0.00483 
0.00485 
0.00486 
0.00487 
0.00487 
0.00487 
0.00488 



0.04412 
0.04344 
0.04317 
0.04278 
0.04265 
0.04248 
0.04240 
0.04232 
0.04225 



0.00424 
0.00422 
0.00423 
0.00423 
0.00423 
0.00424 
0.00425 
0.00426 
0.00426 



04413 
04345 
04318 
04279 
04266 
04249 
04241 
04232 
04226 



00422 
00421 
00421 
00421 
00422 
00423 
00423 
00424 
00424 



an 



0.063, 



ai2 



0.0047 



0.06571 
0.06571 
0.06515 
0.06484 
0.06452 
0.06433 
0.06420 
0.06409 
0.06402 



0.00467 
0.00473 
0.00475 
0.00477 
0.00478 
0.00479 
0.00479 
0.00480 
0.00480 



06572 0, 

06572 0, 

06517 0, 

06486 0, 

06453 0, 

06435 0, 

06422 0, 

06411 0, 

06404 0, 



00467 
00473 
00475 
00477 
00478 
00479 
00480 
00480 
00480 



TABLE II. Phase shifts (in radians) oie+ +H- > e+ + H 
elastic scattering below the positronium formation threshold. 



Ref. Q Ref. Ref. [§§ Ref 



k 



This work 



0.1 


0.1483 


0.152 


0.149 


0.149 


0.1480 


0.2 


0.1877 


0.188 


0.188 


0.189 


0.1876 


0.3 


0.1677 


0.166 


0.166 


0.169 


0.1673 


0.4 


0.1201 


0.118 


0.120 


0.121 


0.1199 


0.5 


0.0624 


0.061 


0.060 


0.062 


0.0625 


0.6 


0.0039 


0.003 




0.003 


0.0038 


0.7 


-0.0512 


-0.053 




-0.050 


-0.0513 
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TABLE III. Partial cross sections (in tvo^) in the 
H{n = 2) - Ps{n = 2) gap (threshold energies 0.7496-8745 
Ry). Numbers 1,2,3 and 4 denote the channels e+ + H{ls), 
e+ + H{2s), e+ + H{2p) and p+ + Ps(ls), respectively. 



Ei{Ry) 



o-ii 



0"12 



0"13 



CTl4 



0.77 Ref. 0.090 0.000702 0.000454 0.00572 

0.77 This work 0.0951 0.000673 0.000331 0.00558 

080 Ref. Hi 0096 0.00115 0.000364 0.00585 

0.80 This work 0.1010 0.00127 0.000371 0.00563 

083 Ref. 13 0.0993 0.00170 0.000885 0.00581 
0.83 This work 0.1063 0.00163 0.000813 0.00566 

084 Ref. 0.101 0.00190 0.00113 0.00580 
0.84 This work 0.1080 0.00173 0.00105 0.00566 




FIG. 1. The short-range part i;*-^-* of the —1/x attractive 
Coulomb potential. 




FIG. 2. The long range part n^'' of the —1/x attractive 
Coulomb potential. 



g_x (E+iZ-z') 



g_Y (z') 



FIG. 3. Analytic structure of g^^ {z — z') gy^ {z') as a func- 
tion of z' with z = E + ie, E < 0, e > 0. The contour C 
encircles the continuous spectrum of hy^ . A part of it, which 
goes on the unphysical Riemann-sheet of gy^ , is drawn by 
broken line. 
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